TRACE FORMULAE FOR PERTURBATIONS OF CLASS S„ 

A.B. ALEKSANDROV AND V.V. PELLER 

Abstract. We obtain general trace formulae in the case of perturbation of self-adjoint 
operators by self-adjoint operators of class S m , where m is a positive integer. In [PSS] 
a trace formula for operator Taylor polynomials was obtained. This formula includes 
6J), the Livshits-Krein trace formula in the case m = 1 and the Koplienko trace formula in 

the case m = 2. We establish most general trace formulae in the case of perturbation 
•^T ' of Schatten-von Neumann class Sm- We also improve the trace formula obtained in 

[PSS] for operator Taylor polynomials and prove it for arbitrary functions in he Besov 
space B™i(R). 

We consider several other special cases of our general trace formulae. In particular, 
we establish a trace formula for mth order operator differences. 
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1. Introduction 



The purpose of this paper is to obtain most general trace formulae for perturbation 
J> ■ of self-adjoint operators by operators of class S m , where m is a positive integer. 

The Livshits-Krein spectral shift function for trace class perturbations of self-adjoint 
operators was introduced by Livshits [L] in a special case and by M.G. Krein [Kr] in 
the general case. It was proved in [Kr] that for a (not necessarily bounded) self-adjoint 



operator A and a trace class self-adjoint operator K there exists a unique function £ in 



O 1 L^M) such that 

o: r 

trace {f(A + K) - f(A)) = / f'{x)i(x) dx (1.1) 

for every function / such that the Fourier transform & f of its derivative belongs to 
L 1 (M). Note that the right-hand side of (1.1) is well defined for every Lipschitz functions 
/. Krein conjectured that for every Lipschitz function / the operator f(A + K) — 
f(A) belongs to the trace class Si and (1.1) holds. It turns out that this is wrong. 
Farforovskaya constructed a counter-example in [F]. 

Later in [Pel] a necessary condition was found. Namely, it was shown in [Pel] that 
if f(A + K) — f(A) € Si for every self-adjoint A and every self-adjoint K in Si, then 
/ locally belongs to the Besov space BL^R) (see § 2 for a brief introduction in Besov 
spaces). This necessary condition also implies that there are Lipschitz functions /, for 
which the condition K £ Si does not imply that f(A + K) — f(A) £ Si. 
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On the other hand, it was shown in [Pel] and [Pe2] that if / belongs to the Besov class 
B^ ol (lR), then the left-hand side of (1.1) belongs to Si and trace formula (1.1) holds. 

We refer the reader to the survey article [BY] for more information about the Livshits- 
Krein trace formula. 

Koplienko considered in [Ko] the case of perturbations by self-adjoint operators of 
Hilbert-Schimidt class S2 ■ With each pair of self-adjoint operators (A, K) such that 
K € S2 he associated a unique nonnegative function 77 in L^R) such that 



, ,. u ,> ( f( A + K )_ f ( A ) - ±f(A + tK) 



4=0 



f"(x)r](x) dx 



1.2) 



for every rational function / bounded on R. 

In [Pe4] the result of Koplienko was improved. It was shown in [Pe4] that if / belongs 
to the Besov space i?^ ol (R), then the operator on the left-hand side of (1.2) belongs to 
S\ and formula (1.2) holds. 

Koplienko also attempted in [Ko] to generalize his results to the case of perturbations 
of Schatten-von Neumann class S m for an arbitrary positive integer m. However, his 
proof for m > 2 was erroneous. 

We also mention here the paper [GPS], in which interesting results related to the 
Koplienko trace formula were obtained. 

In [PSS] it was shown that for every positive integer m and for every pair (A, K) of 
self-adjoint operators with K G S m , there exists a unique function i] m in L X (R) such 
that the following trace formula holds: 



trace [ST^f 



for functions / satisfying the conditions 



f< m \x)r ]m (x)dx 



< j < m, 



;i.3) 



;i-4) 



where the operator Taylor polynomial ^YrJ ls defined by 



^) f ^ f{A + K) - f{A ) 



A,K- 



d 
~dt 



}{A + tK) 



t=o 



(m 



1 d m ~ l 

^1-7 T HA + tK) 



t=o 



The function r/ m is called the spectral shift function of order m. 

Note that earlier partial results were obtained in [DS], [Ski], and [Sk2]. 

In § 6 of this paper we obtain most general trace formulae that include the trace 
formula for operator Taylor polynomials as a special case. 

Moreover, we improve the result of [PSS] for operator Taylor polynomials. We prove 
in § 7 that trace formula (1.3) holds under much less restrictive assumptions on /: we 
show that it holds for all functions / in the Besov space -B™ 1 (R). At the same time 
we establish our general trace formulae also for arbitrary functions in the Besov space 

££i(R). 



In § 5 we establish a formula that expresses the operator Taylor polynomial -^Xk^ 
in terms of a multiple operator integral. This formula will not be used to obtain trace 
formulae. However, we believe that it is of independent interest. 

In § 2 we give a brief introduction to Besov spaces and in §3 we introduce multiple 
operator integrals. Finally, in § 4 we prove two theorems on continuous dependence of 
multiple operator integrals on the corresponding self-adjoint operators that will be used 
to obtain our main results. 

Throughout the paper we use the notation m for Lebesgue measure on R. 



2. Besov classes 

The purpose of this section is to give a brief introduction to the Besov spaces that 
play an important role in problems of perturbation theory. 

Let < p, q < oo and s £ R. The homogeneous Besov class Si (R) of functions 
(or distributions) on R can be defined in the following way. Let w be an infinitely 
differentiable function on R such that 



w > 0, supp w C 



1 

r 2 



and w(x) = 1 — w ( — J for x G [1, 2]. 



We define the functions W n and W n on R by 

&W n {x) = w{£>), &Wl{x) = &W n {-x), neZ, 
where & is the Fourier transform: 

{*f)(t)= I f(x)e- M dx, feL\ 

JR 

With every tempered distribution / 6 J^'(R) we associate a sequences {/ n }nez, 

fn = f*W n + f*Wl (2.1) 

Initially we define the (homogeneous) Besov class S£„ (R) as the set of all / € J^"(R) 
such that 

{2 ns \\f n \\ LP } n& e£i(Z). (2.2) 

According to this definition, the space 5i_(R) contains all polynomials. Moreover, the 
distribution / is defined by the sequence {/ n }nez uniquely up to a polynomial. It is easy 
to see that the series J2 n >o /« converges in ^'(R). However, the series ^2 n<0 f n can 

_ • v^ (r) 

diverge in general. It is easy to prove that the series Xm<o /« converges on uniformly 

R for each nonnegative integer r > s — 1/p if q > 1 and the series ^2 n< Q fn converges 
uniformly, whenever r > s — 1/p if q <1. 

Now we can define the modified (homogeneous) Besov class &L(R). We say that a 

distribution / belongs to B^(R) if {2 ns \\f n \\ LP } n& G F{Z) and f^ = Enez/n m the 
space J?"(R), where r is the minimal nonnegative integer such that r > s — 1/p in the 
case g > 1 and r is the minimal nonnegative integer such that r > s — 1/p in the case 
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q < 1. Now the function / is determined uniquely by the sequence {/ n } n6 z up to a 
polynomial of degree less that r, and a polynomial ip belongs to BpJM) if and only if 
degip < r. 

Besov spaces BpJM.) admit equivalent definitions in terms of finite differences. We 
give such a definition in the case of Besov spaces -B™ 1 (M), with which we mostly deal in 
this paper. 

For t € K, we define the difference operator At by 

{A t f)(x) d ^ f(x + t)-f(x). 
A function / belongs to i?™ 1 (M) if and only if 

|Ar i+1 /||L- , , , |/(a;)| 

1 - • " dt < oo and lim - u \ \{ = 0. 

|i| 1+m |x|-^oo(l + |x|) m 

We refer the reader to [Pee] and [Pe3] for more detailed information on Besov spaces. 



3. Multiple operator integrals 

In this section we give a brief introduction to the theory of multiple operator integrals. 
Double operator integrals appeared in the paper [DK] by Daletskii and S.G. Krein. 
However, the beautiful theory of double operator integrals was developed later by Birman 
and Solomyak in [BS1], [BS2], and [BS4]. 

We are not going to define double operator integrals J J 3>(x, y) dE\ (x) T dE2(y) in the 
case of Hilbert-Schmidt operators T that was the starting point in [BS1], [BS2], and 
[BS4]. We use the approach based on (integral) projective tensor products. In the case 
of bounded or trace class operators T this approach is equivalent to the approach of 
Birman and Solomyak, see [Pel]. Moreover, we start with the definition of the more 
general notion of multiple operator integrals. Multiple operator integrals were defined 
in terms of integral projective tensor products in [Pe5] (see earlier publications [Pa] and 
[St], where multiple operators were defined under much more restrictive assumptions). 

To simplify the notation, we consider the case of triple operator integrals. The def- 
inition for general multiple operator integrals is the same. Let ($?,Ei), (^,£2), and 
(Z, £3) be spaces with spectral measures E\ and E2, and £3 on Hilbert spaces 3%\ and 
J^2, and ^3. Suppose that a function $ on if x f x 2 belongs to the •projective tensor 
product L 00 {E 1 )®L 00 {E 2 )®L 00 {E Z ) of L°°{Ek), L°°(£ 2 ), and L°°(£ 3 ) (i.e., $ admits a 
representation 

®( x i V) = J2 l Pn(x)^n(y)Xn(z), (3.1) 

n>0 

where (p n € L°°(Ei), ip n € L°°(E2), and Xn S L co {E^) are functions such that 

^2 IIVnllzHIVVilliHIXnllL 00 < Oo). (3.2) 

n>0 
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Then for arbitrary bounded linear operators T\ : ,^ 2 — ¥ J%?i and T 2 : J£§ — > ^2, we put 
I f U(x,y,z)dE 1 (x)T 1 dE 2 (y)T 2 dE 3 (z) 




X 9 Z 



ip n dE 2 \T 2 / Xn dE 3 

9 / \Z 

It was shown in [Pe5] (see also [ACDS] for a different proof) that the above definition does 
not depend on the choice of a representation (3.1). For <!> G L°° (Ei)(3L°° (E 2 )(g)L°° (E 3 ) , 
its norm is, by definition, the infimum of the left-hand side of (3.2) over all representations 

(3.1). 

We can enlarge the class of functions <!>, for which multiple operator integrals can be 
defined by considering integral projective tensor products. This approach for multiple 
operator integrals was given in [Pe5]. Again, we consider here the case of triple operator 
integrals; the general case can be treated in the same way. 

We say that a measurable function $ on JT x 'W x Z belongs to the integral projective 
tensor product L°° (Ei)<g>iL°° (E 2 )(3iL°° (E 3 ) if $ admits a representation 



®(x,y,z)= / ( p(.x,uj)ip(y,uj)x(z,ui)d>t(uj), (3.3) 

Jn 

where (f2, xr) is a measure space with a c-finite measure, 99 is a measurable function on 
SC x J), tp is a measurable function on <3f x 0,, \ is a measurable function on Z x Vt, and 

||v(-,w)lk-(Ei)IIV'(-,w)lk-(E 2 )llx(-,w)|| L oo (fi 3 ) dx(u;) < 00. (3.4) 

a 

For a bounded linear operator T\ from M'i to J^i, a bounded linear operator T 2 from 
J#3 to J%, and a function $ in L TO (£i)<g) i L TO (£ 2 )®iI' 00 GE3) of the form (3.3), we put 

f f I 3>{x,y,z)dE l {x)T l dE 2 {y)T 2 dE 3 {z) (3.5) 



X 9 z 




J^(y,uj)dE 2 (y)\T 2 l J 
.9 I \Z 



ip(x,u)dE 1 (x) I Ti I / ip(y,uj)dE 2 (y) ) T 2 \ \ x{z,u)dE 3 (z) } dx{u). 
X / \9 

Again, the above definition does not depend on the choice of a representation (3.3) 
(see [Pe5]). The norm ll^ll^ooA.^oo^.^oo is defined as the infimum of the left-hand side of 
(3.4) over all representations (3.3). 

It is easy to see that the following inequality holds 



I f I 3>(x,y,z)dE 1 (x)T 1 dE 2 (y)T 2 dE 3 (z] 



X 9 z 



< ll $ lli<»«i i I,°°«> i L°°ll r i|l - II^H- 



It is also easy to see that if T\ € S p and T 2 G S 9 , and 1/p + 1/g < 1, then the triple 
operator integral (3.5) belongs to S r and 



J J [$(x,y, z) dE x {x)T x dE 2 (y)T 2 dE 3 {z] 



X 9 Z 



^ \\^\\L°°®iL' x ®iL ca 11-^1 lISp ' ||-*2||Sgj 



where 1/r = 1/p + 1/q. 

Note that in [JTT] Haagerup tensor products were used to define multiple operator 
integrals. However, it is not clear whether this can lead to stronger results in perturbation 
theory. 

For a function / on R the divided differences T) k f of order k are defined inductively 
as follows: 

2)°/ = /; 

if k > 1, then in the case when ti, t 2 , • ■ • , tfc+i are distinct points of R, 

/^fc ,w n dcf (S fc-1 /)(^i, • • • , a;fc-i, Xfc) - (3? fc_1 /)(^i ! ■ ■ ■ , x fe _i, Xfc+i) 
(i) /J(»i,--- ,X k+ i) = 

x k ~ x k+l 

(the definition does not depend on the order of the variables) . We put 

If / £ C fe (R), then 2) fc / extends by continuity to a function defined for all points 
x 1 ,x 2 ,--- ,X k +l- 

It can be shown that 

m+l fc— 1 m+1 

(2) m /)(xi,...,x m+ i) = 5^/(ar fc ) n^* - ^')" 1 II (^fc-^')" 1 - 

fc = l j = l jr'=£+l 

It follows from the results of Birman and Solomyak [BS4] that if A is a self-adjoint 
operator (not necessarily bounded), K is a bounded self-adjoint operator, and / is a 
continuously differentiable function on R such that 2)/ € L 00 (£' J 4_|_x) ( X'i-^ 00 (-EA), then 

/(A + if) - /(A) = J! (2)/) (x, y) dE A+K (x)K dE A (y) (3.6) 

MxR 

(throughout this paper Ea stands for the spectral measure of A). 

It was shown in [Pe5] that if / is a bounded function on R such that its Fourier 
transform &f is supported on [— a, a], then then 

$ro/ G L 00 (R)®i---<8>iL 00 (R) (3.7) 

V v ' 

m+l 

and 

P^Wl-v^l- < const a m \\f\\ L ~. (3.8) 

This implies (see [Pe5]) that if / belongs to the Besov space S™ 1 (R), then (3.7) holds 
and 

ll^/IU^L- < const ||/||b SiW . (3.9) 



It was also proved in [Pe5] that if / £ -B^ 1 (R), A is a self-adjoint operator and K is 
a bounded self-adjoint operator, then the function 

t*->f(A + tK) 

has m derivatives in the operator norm and 

d" 1 



itm ttA + tK) 



t=0 

m+1 



J---j{® m f){xi,--- ,x m+1 )dE A ( Xl )K ■ ■ ■ KdE A (x m+1 ). (3.10) 



Strictly speaking under the condition / € i?™ 1 (M) lower derivatives 

d k 
—f(A + tK), k<m, 

do not necessarily exist. To get the mth derivative we should define the left-hand side 
of (3.10) by 

— f{A + tK) = J2^fn(A + tK), 

neZ 

where the functions f n , n £ Z, are defined by (2.1). 



4. Continuous dependence of multiple operator integrals 

We have already seen (see (3.10)) that multiple operator integrals that involve divided 
differences of an arbitrary order play and important role in perturbation theory. In this 
paper we are going to consider other multiple operator integrals that involve divided 
differences. 

The purpose of this section is to show that such multiple operator integrals depend 
continuously on the corresponding self-adjoint operators. We establish two similar re- 
sults. The first one establishes continuous dependence in the operator (or trace class) 
norm, while the second one deals with continuous dependence in the strong operator 
topology. 

We are going to use a representation of T) m f as an element of the integral projective 
tensor product of L°° spaces. This representation was obtained in [Pe2] for m = 1 and 
in [Pe5] for m > 2. To simplify the notation, we formulate the result for m = 1 and 
m = 2. Similar formulae hold for m > 2. 
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Let / be a bounded function on M. such that the support of its Fourier transform &f 
is a compact subset of (0, oo). Then (see [Pe2]) 



VDf)(x,y) = ± jj (^f)(s + t)e isx e [t ydsdt 

R+xR+ 



If — {^f)(s + t)e isx e ity dsdt 



2vr J J s + t 

.xR+ 



+— ff -^(^f)(s + t)e iSX e lty dsdt 
2vr J J s + t 



.XJ1 



{(S- t f)*q t )(xy t «dt + i / e isx ((S- s f)*q s )(y)ds 



where (S t f)(x) = e ltx f(x) and qt is the distribution defined by 

(^qt)(s) = max { 1 - r -^— , 

One can prove that for t > 0, qt is a (complex) measure whose total variation does not 
depend on t, see [Pe2]. 
Hence, 

(2>f)(x,y) = i [ ((S- t f)*q t )(x)e i »dt + i f e isx {(S~ s f) * q s )(y) ds (4.1) 

for every / € L°°(IR) such that the support of its Fourier transform J^~/ is a compact 
subset of (0, oo). Moreover, if supp J^/ C [0, a], a > 0, then 

/ \\{(S-'f)*q a )\\ LOB dt< const (7||/|| L =o. (4.2) 

In the same way we have 



' m+l 



:m r r I ^ \ 

(V m f)(x 1 ,x 2 ,...,x m+1 ) = —j---j(^f) [J^tj e^V^.-.e^^cft 



w ! 

-1 



m+l 

£ 

i =1 o o U= ] / U =1 



: ~ IZ / " " " / *J 5Z *J W) 2 *J e itlXl e it2;C2 . . . e iim+ia;m + 1 eft 



for every function / £ C(M) such that &f € L 1 (M) and suppJ^/ C [0, oo). Applying 
this formula for m = 2, we obtain (see [Pe5]) 

(&f)(x,y,z)=- J J {(S- t - u f)*q t+u )(x)e it ye iuz dtdu 

R+xR+ 



II e isx ((S- S ~ u f) * q s+u ) (y)e iuz ds du 



e isx ^y ((«S" S -*/) * q a +t) (*) ds dt (4.3) 



for every / € L°°(K) such that the support of its Fourier transform &f is a compact 
subset of (0, oo). 

Moreover, if supp J?/ C [0, a], a > 0, then 

| ((£-'-"/) *qt+u)\\ Loo dtdu < consta 2 ||/||z~. (4.4) 



Theorem 4.1. Let m be a positive integer and let f € 5J2 1 (R). Suppose that 
Ki, K 2 ,--- , K m are bounded linear operators, Ai, A 2 ,--- ,A rn +i are (not necessarily 
bounded) self-adjoint operators, and {.Aijj , {^4 2;J } > n ) ' " ■> {Am+i,j} >n are sequences 
of self-adjoint operators such that 



lim || A\ - AijH = lim ||A 2 - A 2 j|| = • • • = lim ||-A m+ i - A m+ i ;J -|| = 0. 

j—>oo " j—toc ' j—^OO 



Then 



lim / ••• / (D m /)(x ir -- ,x m+1 ) dE A (x 1 )K 1 dE A2 Ax 2 )K 2 ■ ■ ■ K m dE A (x m+1 ] 

— >oo J J >J .J i .J 



m+1 



= /"■■■ |(D m /)(x 1 , • • • , x m+1 ) d£ Al (xi)Ki dE M {x 2 )K 2 ■■■K m dE Am+1 (x m+1 ) (4.5) 

m+l 

in the operator norm. 

Moreover, if K±, K 2 , ■ ■ ■ ,K m G S m , then the limit in (4.5) exists in the norm, of S\. 

Proof. The proof is similar to the proof of Theorem 5.1 in [Pe5]. For simplicity, we 
give the proof in the case m = 2. For arbitrary m, the proof is the same. 

As usual, it suffices to prove the result in the case when supp &f is a compact subset 
of (0, oo). We are going to use formula (4.3). Put 

$(x 1 ,x 2 ,x 3 ) = // {(S- t - u f)*q t+u )(x 1 )e itX2 e iuX3 dtdu. 



We have to show that 



Ihn /// $(x 1 ,X2,x 3 )dE Al .(x 1 )KidE A2 .(x 2 )K 2 dE A3 .(x 3 ) 

$(xi,X2,x 3 )dE Al (xi)KidE A2 (x2)K 2 dE A3 (x 3 ). 



The proof for the other two terms in (4.3) is exactly the same. 
Clearly, 



$(zi, x 2 , x 3 ) dE Alj {x{)Ki dE A2j (x 2 )K 2 dE A . Aj (x 3 ) 
jj {{S^ u f) * q t+u ){A lt3 )K 1 e tA ^K 2 e uA ^ dtdu 



and 



$(ari, x 2 , x 3 ) dE Al (x 1 )K 1 dE M (x 2 )i^2 dE Az (x 3 ) 
J J ((5-*- u /) * q t+u )(A 1 )K 1 e itA "K2e iuA " dtdu. 



It is easy to see that it suffices to show that 

lim ||((5-*-«/) * qt+u)(A^) - {(S' l - U f) * gt+uXAOll = 0, 



lim ||e iiA ^ - e iM2 || = lim ||e iMA ^ - e iuA3 || = 

j—too j— >°° 

for every t, u > 0. However, this is obvious, since the functions (<S~*~ M /) * qt+ u and 
x i— ?• e ttx are operator Lipschitz. 

The same reasoning shows that if K±, K 2 , • • • , K m € S m , then the limit in (4.5) exists 
in the norm of Si. I 

Theorem 4.2. Lei m be a positive integer and let f € i?™ 1 (M). Suppose that 
K\, K2, ■ ■ ■ , K m are bounded linear operators, A\, A2, ■ ■ ■ ,A m+ \ are (not necessarily 
bounded) self-adjoint operators, and {^4ij} >0 , {^2,j} >0 , ■ ■ ■ , {An+i,j} > are sequences 
of bounded self-adjoint operators such that 



lim \\AkX — Ak i'x|| = 
10 






for k = 1,2,... , m + 1 and x m £/te domain of A/,. Then 
lim / • • • / (® m /)(^i, ■ ■ ■ , arm+i) d^u (^i)-^i ^A 2J (^2)^2 ■ ■ ■ K m dE Am+1}] (x m +i) 

JTl+l 

(D m /)(xi, . . . , x m+ i) d^i (a;i)ifi dE A2 (x 2 )K 2 . . . K m dE Am+1 (x m+1 ) 

m+l 

in the strong operator topology. 

We need two lemmata. The one is Lemma 8.4 in [API]. 

Lemma 4.3. Let f be a bounded continuous function on R. Suppose that A is a 
self-adjoint operator (not necessarily bounded) and {Aj}j>q is a sequence of bounded 
self-adjoint operators such that 



Then 



lim \\AjU — Au\\ = for every u in the domain of A. 



lim f(Aj) = f(A) in the strong operator topology. 

j->oo 



The second lemma is an obvious operator-valued version of the Lebesgue dominated 
convergence theorem. 

Lemma 4.4. Let (fi, /i) be a measure space. Let {&j}j>o be a sequence of weakly 
measurable operator-valued functions. Suppose that sup ? - ||$j(w)|| £ L 1 ^) and $(w) = 
lim $j(uj) in the strong operator topology for /i-almost all u. Then 

lim / ®jd[i= / $ dfi 

in the strong operator topology. 

Proof of Theorem 4.2. Clearly, it suffices to prove the result in the case when 
supp^/ is a compact subset of (0, 00). By Lemma 4.3, 

lim exp(itA kJ ) = exp(fL4 fc ) and lim ((<S - */) * Qt){A Kj ) = ((<S - */) * qt)(A k ) 

j— >oo j— >oo 

in the strong operator topology for all t > 0. Now the result follows from Lemma 4.4 
and the integral representation for T) m f as an element of projective tensor product of 
L°° spaces, see the discussion in the beginning of this section. ■ 



5. A formula for operator Taylor polynomials 
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In this section we obtain a formula for operator Taylor polynomials in terms multiple 
operator integrals. We are going to prove that for / € i?™ 1 (M), the following formula 
holds 



[■■■f(T) m f)(x 1 ,--- ,x m+1 ) dE A+K ( Xl )K dE A (x 2 )K ■ ■ ■ K dE A (x m+1 ). (5.1) 



ar( m ) f 
' y A,KJ 

ra+l 

We believe that formula (5.1) is of independent interest though it will not be used to 
prove the the main results of the paper. 

Here A is a (not necessarily bounded) self-adjoint operator and K is a bounded self- 
adjoint operator. 

Recall that 

,(m) p def 



nilf = f^ + K)-f{A) 



d 



dr 



TO— 1 



—f{A + tK) - — -f(A + tK) . (5.2) 



t=o (m — 1)! dt 



t=o 



Actually, it is not true that each term of the right-hand side of (5.2) exists for functions 
/ in i?™ 1 (M). However, we explain in this section how to define the right-hand side of 
(5.2) for all functions / in B^iW- 

First, we establish (5.1) for bounded functions / whose Fourier transform has compact 
support in (0, oo). 

Theorem 5.1. Let f be a function in L°°(R) such that the support of ^ f is a compact 
subset of (0, oo). Then (5.1) holds. 

To prove Theorem 5.1, we need the following result. 

Lemma 5.2. Let f be a function in L°°(M) such that the support of & f is a compact 
subset of (0, oo). Then for every m > 2, 

y. • J (ST" 1 /) (m, ■■■ ,x m ) dE A+K ( Xl )KdE A {x 2 )K ■ ■ -KdE A {x m ) 

m 

~ / ' / ( S)Tn ~ 1/ ) (a;i ' ' ' ' ' Xm) dE A {x l )KdE A (x 2 )K-- -KdE A (x m ) 
:i 
= [■■■[ (® m f)(x!, • • • , x m+1 ) dE A+K (x 1 )KdE A (x 2 )K ■ ■ -KdE A (x m+1 ). 

m+l 



Proof. To simplify the notation, we give the proof in the case m = 2. The proof in 
the general case is exactly the same. 
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Consider first the case when A is bounded. We have 




jj{T>f)(x,y)dE A+K (x)KdE A {y) - jj{T>f)(x,y)dE A {x)KdE A (y) 
[[(Vf){x,z)dE A+K (x)KdE A {z) - JJ(Vf)(y,z)dE A (y)KdE A (z) 
(£/)(x, z) dE A+K (x) dE A (y)K dE A (z) 
(2)/)(y, z) dE A+K {x) dE A {y)KdE A {z) 

III (( I)2/ )^' y ' z ^> ( x ~ y ) dE A+K{x) dE A (y)KdE A {z) 
£> 2 /) (x, y, z)) x d£? A +A-(x) d£ A (y)K dE A (z) 
£> 2 /) (x, y, *)) y dE A+K {x) dE A {y)K dE A {z) 
If {V 2 f)(x,y,z)dE A+K (x)(A + K)dE A (y)KdE A (z) 

~ III (® 2 f)( x >y>^ dEA + K W AdEA ^ KdEA W 

HI {® 2 f)(x,y,z)dE A+K (x)K dE A (y)K dE A (z). 




In the case of unbounded A, we put 
Cieariy, 



Aj^AEaH-JJ]). 



lim ||AjU — Au|| = for every u in the domain of A. 
Since each operator Aj is bounded, we have 

JJ(J)f)(x,y)dE A . +K (x)KdE A .(y) - J J \®f)(x,y) dE A] {x)K dE A] {y) 

(D 2 /) (x, y, z) dE Aj+K (x)K dE Aj {y)K dE Aj (z). 

The result follows now from Theorem 4.2. ■ 

Proof of Theorem 5.1. We proceed by induction. It was proved in [Pel] that 
identity (5.1) holds for m = 1 and for functions / in the Besov space B^ ol (R), and so it 
holds for functions / satisfying the hypotheses of Theorem 5.1. 
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To pass from m — 1 to m, we use Lemma 5.2 and the formula 

-./(-1+/./V) 



„«„,-/(•"+«) 



=m! !■ ■ j (» m /)(n, ■ ■ ■ , x m+ i) dE J+K (x,)A'6!Ej(x 2 )A' ■ ■ ■ A'<iE 4 (i- m+1 ). 

771+1 

This formula was proved in [Pe5] for / G S™ 1 (]R), and so holds for functions satisfying 
the hypotheses of Theorem 5.1. I 

We can extend now formula (5.1) to the case of arbitrary functions / in i?™ 1 (M). 
As we have already mentioned, for / £ S™ 1 (IR), each term of the Taylor polynomial 

^Xxf is not necessarily defined. However, we can define ^Xxf ^ ro / e -^ooi(^) by the 
following formula: 

oo 
^(m) , def ^ ^(m) f ,_ „\ 

' y A,KJ — 2-*i ' y A,KJn, KP-Z) 

n= — oo 

where the functions f n are given by (2.1). Note that it follows from (3.8) that the right- 
hand side of (5.3) converges absolutely for all functions / in 5^ 1 (R). It is also easy to 
see that the right-hand side of (5.3) does not depend on the choice of the function w in 
the definition of Besov spaces, see § 2. 

Theorem 5.3. Let f G B™^). Then (5.1) holds, where &j$f is defined by (5.3). 

Proof. By Theorem 5.1, (5.1) holds for functions / in L°°(R) such that suppJ^/ 
is a compact subset of (0, oo). In the same way it can be proved that (5.1) holds for 
functions / in L°°(R) such that supp^"/ is a compact subset of (— oo, 0). Thus (5.1) 
holds for each function f n . The result follows now from Theorem 5.5 of [Pe5]. I 



6. The general result 

In this section we establish most general spectral formulae for perturbations of class 
S m (see Theorem 6.5 below). We show in the next section that the trace formula for 
operator Taylor polynomials is a special case of Theorem 6.5. 

Let A be a self-adjoint operator and let K be a bounded self-adjoint operator. Put 

dcf 

At = A + tK for i£K. Let / be a bounded continuous function on R. and let m be a 
positive integer m. We can consider the following finite differences: 

TO , N. 

(A%f)(A)^Y,(-ir-i I™) f(A + j K). 

3=0 V J 

It turns out that the finite differences (A^/) (A) can be defined for functions / in B™-± 
by the formula 

(a£/)(V=E( a M(^ 
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where the functions f n are defined by (2.1). Moreover, it was shown in [API], Lemma 
4.3 that 

(A?/) (A) 

=m! /■ • ■j{T) m f){x l , ■■■ , x m+1 ) dE A { Xl )K dE A+K {x 2 )K ■ ■ ■ K dE A+mK {x m+1 ). (6.1) 

m+l 

Strictly speaking, formula (6.1) was proved in [API] for bounded self-adjoint operators 
A. However, it is easy to see that the approximation procedure used in the proof of 
Lemma 5.2 in this paper also works to extend formula (6.1) to the case of unbounded A. 
Recall that it was shown in [Pe5] that for every function / in S™ 1 (IR), the function 
t h- > f(At) has mth derivative in the operator norm if we define it by 

d m h f d m 



and 



dp 



-ml 



-f(A t 



m+l 



t=s 



(D m /)(xi,--- ,x m+1 )dE As {x 1 )KdE Aa (x 2 )K---KdE As (x m+1 ), (6.2) 



where E s is the spectral measure of A s . 

The following estimate was obtained in [PSS] for functions / satisfying (1.4). We 
extend it to the class B™ ± ( 



Theorem 6.1. Let f £ £™i(K) and K G S m . Then 

d m 



trace 



dV 



;f(At 



L x 



< const ||/ (m) |L 

— II J IIL C 



\K\ 



(6.3) 



Proof. Clearly, it suffices to show that 

d 



trace 



dV 



-f(A 



t=o 



< const ||/ (Tn) || r J|K||s 

— IK II L " " °rr 



Let us first prove the result in the case when supp J^~/ is a compact subset of (0, oo). 
Let $ be a function in C°°(R) such that 3>(0) = 1 and ^<& is a nonnegative infinitely 

differentiable function with a compact support. For e > 0, we put f £ {x) = &(ex)f(x). 
Then supp &f £ is a compact and 

&f e £L\WL)nC°°(WL). 

Hence, by Theorem 2.1 of [PSS], 



d m 
traCe ( dt^ fe(At: 



t=o 



< const ||/J m ) \\ TOO \\K\ 



15 



It is easy to see that 



I f( m )|l < n\\ f( m )\\ 

\J£ II L°° 1 1*' IlL 00 ' 



where C depends only on $ and m. Moreover, supp # f e is a compact support of (0, oo) 
for sufficiently small e and lim £ ^o /e = / m the space i?™ 1 (lR). It follows from (6.2) and 
from (3.8) that 



d" 



-(f-fe)(At) 



Hence, 



dt< 



d m 
lim trace I — — f e (Ai 

e^O \dt m V ' 



< const 



h 



eWB 71 



\K\ 



Si 



t=0 



d m 
trace ( - — f (At 



t=o 



(6.4) 



han ' ' dt^^ At/ 



< const 1 1 / 



(m) I 



\K\ 



If suppJ^/ is a compact support of (— oo,0), the proof of (6.4) is the same. 
It follows from (6.4) that for / in B^iR), 

d n 

' - ./„(.•!/) 

t=0 / 

where the functions f n are defined by (2.1). This implies (6.3). ■ 

Theorem 6.2. Let m be a positive integer and let f € S™ 1 (M). Then 
d m 



dV 



-f(At 



= limh- m (A m 

t=s /i-s-0 



hK f)(As). 



(6.5) 



Proof. The result follows from (6.2), (6.1) and the equality 

lim / ••• /(D m /)(xi,--- ,x m+1 ) dE A {x 1 )K dE A+hK (x 2 )K ■ ■ ■ K dE A+mhK (x m+1 ) 



m+l 



= j ■■ y (D m /)(x 1; • • • , xm+i) dE A { Xl )K dE A (x 2 )K ■ ■ ■KdE A (x m+1 ) 

m+l 

which is a special case of Theorem 4.1. ■ 
Consider now the case when K £ S m . 

Theorem 6.3. Let f G B^R) and K € S m . Then the limit 

Kmh-™(Af K f)(A s ) 

exists in the norm of Si for every s. Moreover, the function 

d m 



S H-> 



d,f 
is a continuous Si-valued function and 

d n 



;f(A t 



t=s 



dt< 



-f(At 



< const 



L°°(Si) 






(6.6) 
(6.7) 
(6.8) 



k> 



Proof. The fact that the limit (6.6) exists in the norm of S\ is an immediate con- 
sequence of Theorem 4.1. The continuity of the function (6.7) also immediately follows 
from Theorem 4.1. Finally, inequality (6.8) follows from (6.2) and inequality (3.9). ■ 

Let Co(R) denote the space of all functions / G C(R) such that limui^^ /(£) = 0. 
Denote by ^#(R) the space of all finite Borel measures on R. We identify in a natural 
way the space ^#(R) with the space (Cq(R))*. 



Theorem 6.4. Let m be a positive interger and let A and K be self-adjoint operators 
such that K G S m . Then for each £ £ R there exists a unique measure v t G ^(R) such 
that 

(im \ r 

— f(A t )\ = / f im) du t (6.9) 

for every f G S™ X (R). Moreover, the map t t-^- v t is a continuous map from 
equipped with the weak-* topology cr(^#(R), Cq(R)). 



£o ^(R) 



del' 



Proof. Note that the set X =* ^^(R) n C (R) is dense in C (R). Theorem 6.1 
implies that there exists a unique measure v% G ^#(R) such that (6.9) holds for every 
/ G X. Theorem 6.3 implies that the mapping £ i— > vt is continuous. 

It remains to establish the equality in (6.9) for all / G i?™ 1 (R). Clearly, it suffices to 
verify that this equality holds on a dense subset of i?™ x (R). Let / be an entire function 

of exponential type a that is bounded on R. Put tp(t) = t~ l sin£. It is easy to see that 



1*1 



lim (</?(££)/(£)) 



(ra) 



, ( m ) 



0, 



iH/ 



lim(^( £ £)/(£))^ = (/(£))^(£), 



and 



sup 
e.t 



{V(et)f(t)) 



(m) 



< CO. 



To complete the proof, it suffices to verify that 



lim trace 

£^0 



df 



■(ip(A £t )f(A t )) = trace 



We have 



d m 

— ( v {A Ei )l(A,)) 



j=0 



df 



d m ~i 



;f(A t 



— 

dP^ et ' dl' n -i 



f(At 



Inequality (5.3) in [AP2] implies the following inequalities: 
< C m £ J ||.K"||'Q and 

S rn 



d3 ,a ^ 



dt m -i 



f(At 



< C rr m ~J\\K\\ m ~ 3 



Hence, 



lim 

£-5>0 



dP 



dP 



M A 



Et) 



fjpn-j 



dt m -3 



jf(At) 



0, J>1. 
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It remains to prove that 

/ d m \ ( d m 

lim trace \^p(A et ) —f(A t )j = trace \^f( A t) I 

which is clear, because lim ip(A £ t) = I in the strong operator topology. ■ 

For pu € ^#(M), we denote by jjS" 1 ' the m-th derivative of /i in the sense of distributions. 
Thus we put 

\ a r [' d m 

f(A t ),^)^(-irj^— f(A t )dn for /eC, 



Let 



def / -. n 
V{A,K} = (-1) 



I \dtJt(t), (6.10) 

Jm. 



the integral exists in ^(R) in the topology cr(^#(lR), Co(R)), because the function t^rvt 
is weak star continuous. 

We prove in the next section that if \i is an absolutely continuous measure (with 
respect to Lebesgue measure), then H{a,k} must also be absolutely continuous. 

The following result is a most general trace formula for perturbations of class S m . 

Theorem 6.5. Let m be a positive integer and let A and K be self-adjoint operators 
such that K £ S m . Suppose that [i £ „#(R). Then for every f £ S™ 1 (R) ; 



trace (f(A t ),^ m A = / /<*") d» {AtK} , (6.11) 

x ' JR 

where ^sa,k\ * s defined by (6.10). 
Proof. We have 

trace (f(A t ),^) = (-l)™jf trace (|^/(A)) ^ 

= ("l) m [ ( [ / M dvt) dp= I f^ dp {AK} . U 

JR \JR / JR 

7. Trace formulae for operator Taylor polynomials and other special cases 

In this section we show that trace formula (1.3) is a special case of our general formula 
(6.11) which allows us to improve the main results of [PSS] and extend trace formula 
(1.3) to the case of functions in the Besov space B'^ )1 (R). We also consider several other 
interesting special cases of formula (6.11). 

Theorem 7.1. Let m be a positive integer and let A and K be self-adjoint operators 
such that K € S m . Consider the the absolutely continuous measure fi defined by 

d ^) = T^iyyt 1 " *) ra " 1 x[o,i](*)*- (7- 1 ) 
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then the formula 

trace (VjJ/) = f f {m) (x)d^ {AK} (x) 
holds for every function f in the Besov space i?™ 1 (IR). 
Proof. It is easy to see that 

where 5 a is the unit point mass at a. Indeed, by Taylor's formula, 

, (1) J£^ + /'<i^> (t)(ft 

PJ j! Jo (m-1)! 

for every <p G C m (lR). It follows that 

^j$/=(/(A),M (m) )- 

The result follows now from Theorem 6.5. ■ 

Corollary 7.2. Under the hypotheses of Theorem 7.1, t/ie spectral shift function rj m 
of order m is the Radon-Nykodym derivative of H{a,k} with respect to Lebesgue measure, 

d^{A,K} _ 
dm 
and trace formula (1.3) holds for every f G 5™ 1 (1R). 

Proof. It follows from Theorem 7.1 and from trace formula (1.3) that 

f (m) (x)^ { A,K}= I f {m \x) Vm { X )dx 

JR 

for every infinitely smooth functions / with compact support. This implies the result. 
■ 

Now we are in a position to prove that H{a,k\ is absolutely continuous whenever jj, is. 

Theorem 7.3. Suppose that jjl is an absolutely continuous measure (with respect to 
Lebesgue measure). Then /j-sa,k} i s a ^ so absolutely continuous. 

Proof. For g G ^(m), we use the notation 

9{A,K} = H{A,K}, 

where /j, is the absolutely continuous measure defined by dfj, = g dm. 

Denote by Jz? the set of all g G L l {m) such that the measure g{A,K} 1S absolutely 
continuous for all self-adjoint operators A and K with K G S m . Clearly, «5f is a closed 
translation and dilation invariant subspace of ^(m). By Theorem 7.1 and Corollary 
7.2, the function i4(l- £) m_1 X[o,i](*) belongs to Jgf. Since /^(l - t) m ~ l dt / 0, it 
follows that J*f = L l (m). ■ 

The following results are interesting special cases of our general trace formula (6.11). 
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Theorem 7.4. Let m be a positive integer. Suppose that A is a self-adjoint operator 
and K is a self-adjoint operator of class S m . Then there exists a function x m € L 1 (M) 
such that 

trace (A%f){A) = [ f^ m \x)K m {x) dx 

for every f € 5^ 1 (R). 

We need the following lemma. 

Lemma 7.5. Let m be a positive integer and let X be a measure on a finite interval 
[a,b]. Suppose that X is orthogonal to the polynomials of degree less than m. Then 
there exists a function g € L°°(R) such that g( m > = X in the sense of distributions, 
supp g C [a, b] and \\g\lL 1 < (m!) _1 (6 — a) m ||A| 



Proof. Without loss of generality we may assume that a = 0. Put 

9 (x)= r {x ~ t)m ~' dx(t). 

Jo (m-1)! 
It remains to observe that g{x) = for x [0, b] and 

m— 1 

1^)1 - ( m -i)! l|A|| -^W 
for x £ (0,6). ■ 

Proof of Theorem 7.4. Put 

A^£(-ir-'( m k, 

By Lemma 7.5, there exists a function g € L°°(R) such that supp g C [0, m] and [i = g( m ' 
It is easy to see that 

(A%f)(A)= [ f(A t )dX(t). 
Hence, by Theorem 6.5 we have 

trace ((A?/) (A)) = / f™ dn {A , K} , 

JR 

where d/j,(t) = g(t) dt. Moreover, the measure H{a,k\ is absolutely continuous, see § 6. 
■ 

Now we state the following generalization of Theorem 7.4. 

Theorem 7.6. Let m be a positive integer and let A and K be self-adjoint operators 
such that K G S m . Suppose that the measure ^7=0 ^j^tj * s orthogonal to the polynomials 
of degree less than m. Then there exists a function ip in L 1 (R) such that 

N 



J=0 

for every f eB^fl). 



trace £ X 3 f{A t] ) = jf / (m) (*Mx) dx 
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The proof of Theorem 7.6 is similar to that of Theorem 7.4 and we omit it. 
Let us now generalize Theorem 7.6. 

Theorem 7.7. Let m and k be nonnegative integers such that < k < m. Suppose 
that the measure J2j=o x j$tj * s orthogonal to the polynomials of degree less than m — k. 
Let A be a self-adjoint operator and K be a self-adjoint operator of class S m . Then there 
exists a function ip G L 1 (R) such that 




d m ° 
trace I \ "A,-- f(A 



H, 



f( m \x)ip(x) dx 



for every f G i?™ 1 (R). 

Proof. By Lemma 7.5, there exists a function g G L°°(R) with a compact support 
such that and g( m ~ m o) = Y.f=o X J 6 t y Then g^ = Y2f=o X J d t™ o) ■ B Y Theorem 6.5, we 



'■j 
have 

N 



ffmo r 



3=0 

where Hia,k} 1S defined by (6.10) with dfj,(t) = g{t)dt. Clearly, the measure ^sa,k} ls 
absolutely continuous with respect to Lebesgue measure. ■ 

We have considered several special cases our general trace formula (6.5). In all those 
cases the measure Hsa,k\ 1S absolutely continuous with respect to Lebesgue measure. 
The following example shows that this is not always the case. 

Example. Suppose that A and K are commuting self-adjoint operators with discrete 
spectra, i.e., A = J2j a j('i e j) e j an< ^ K = Ylj X j('' e j) e j' where {ej}j is an orthonormal 
basis and the ctj and \ij are real numbers. By Theorem 6.4, 



*» (£/<*>) -j^-'* 



for a unique measure v% on R. Note that this is a special case of trace formula (6.5). It 
is easy to verify that 

MZfM) = ^ taCe (/(^)) = £ X ?f (m \«3 + *A,-) 

3 

for every / G C m (R) with /( m ) G L°°(R). Thus 

Ut = J2 X T 6a j+ tx j- 
j 

Hence, v t is a discrete measure. 

This example shows that in trace formula (6.11) the measure H{a,k\ does not have to 
be absolutely continuous in general. 
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